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ABSTRACT 

In  this  paper,  modified  strength  of  materials  (MSM)  and  strain  energy-based  approach  have  been  used  to  determine 
effective  electro-mechanical  properties  of  piezoelectric  fiber  reinforced  composites  (PFRC).  Based  upon  these  methods, 
micromechanics  model  has  been  formulated  using  a  rectangular  representative  volume  element  (RVE)  which  contains 
both  fiber  as  well  as  matrix  phases.  Result  obtained  has  been  compared  with  strength  of  materials  (SM)  approach  to  the 
similar  problems  available  in  literatures.  Results  suggest  that  there  is  wide  variation  in  estimation  of  properties  obtained 
through  strength  of  materials  method  to  the  methods  adopted  in  this  work.  Strength  of  materials  (SM)  method  in  most  of 
cases  overestimates  or  underestimates  the  actual  result.  Though  each  of  such  approaches  has  its  own  limitations,  the 
adopted  method  in  this  paper  refines  result  obtained  through  strength  of  materials  (SM),  as  it  takes  into  account  of  the 
several  conditions,  which  has  not  been  dealt  while  analysing  with  strength  of  materials  method.  Both  modified  strength  of 
materials  and  energy  approach  reduces  aberration  generated  in  result,  due  to  generation  of  unequal  strain  in  matrix  and 
fiber  phases  as  well  as  Poisson ’s  ratio  mismatch  between  fiber  and  matrix  phases. 
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Nomenclature 


Symbol 

af 

af  ( i  =  1  to  4) 
a™  (i  =  1  to  4) 
At  (t  =  1  to  3) 
Bt  (i  =  1  to  7) 
C 
D 
e 
E 

Ft  (i=l  to  4) 
Pi  (i  =  1  to  18) 
G;  (t  =  1  to  4) 
i,j,  k,l,m,n 
/;  (t  =  1  tO  4) 
Ji  (i  =  1  to  4) 

RiJ 

a 

e 

K 


Description 

Area  of  cross-section  of  original  fiber 

Constants 

Constants 

Constants 

Constants 

Stiffness 

Dielectric  displacement 
Piezoelectric  stress  coefficient 
Electric  field 
Constants 

Experimental  Constants 
Constants 

Indices  representing  field  quantities 

Constants 

Constants 

Ratio  of  piezoelectric  stress  coefficient  of  PFRC  to  piezoelectric  fiber 

Stress 

Strain 

Dielectric  Constant 


Unit 


N/nr 

C/m2 

C/m2 

V/m 


N/nr 

C/Vm 
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1.  INTRODUCTION 

Piezoelectric  materials  have  an  excellent  property  of  doing  reciprocity  of  energy  between  the  mechanical  domain  to 
electrical  domain  and  vice  versa.  This  makes  piezoelectric  ceramics  such  as  PZT  (Lead  Zirconium  Titanate)  very  attractive 
materials  towards  sensors  and  actuators  applications.  Piezoelectric  materials  have  numerous  applications  in  conventional  as 
well  as  advanced  engineering  applications  like  aerospace,  microphones,  ultrasonic  transducers,  automatic  control  of  remote 
explorations.  These  wide  applications  make  piezoelectric  materials  a  class  of  materials  that  requires  attention  for  research 
and  development  toward  miniaturisation  of  properties  of  these  ceramics.  While  using  in  pure  form,  these  materials  have 
certain  limitations,  like  low  piezoelectric  constants,  shape  control  and  high  specific  acoustic  impedance,  this  can  be 
overcome  by  combining  these  materials  with  non-piezoelectric  materials.  This  results  into  piezoelectric  composites  bearing 
superior  mechanical  and  electrical  properties  to  the  previous  one. 

The  study  of  electromechanical  properties  of  piezoelectric  fibre  reinforced  composites  (PFRC)  has  drawn  a  great 
deal  of  attention  in  recent  years  in  form  of  mathematical  modelling  of  behaviour  of  these  composites.  Various  numerical 
methods  as  well  as  analytical  methods  has  been  used  to  describe  the  behaviour  of  such  elements,  though  these  methods 
have  limitations  of  their  own.  Currently,  there  has  been  increasing  research  interest  in  developing  analytical  models  to 
determine  the  effective  electro-mechanical  properties  of  such  class  of  materials.  In  micromechanical  approach,  the  overall 
properties  of  the  composite  are  said  to  be  volume  average  sum  of  the  properties  of  a  unit  representative  volume  element 
(RVE)  constituted  of  both  fiber  and  matrix  surrounding  it.  Effective  properties  can  be  derived  from  the  results  obtained  by 
the  analysis  of  the  RVE.  A  number  of  such  efforts  for  determining  effective  electro-mechanical  properties  of  such 
composites  [6-8]  have  been  reported  in  the  literature. 

All  of  such  models  have  many  assumptions,  which  are  common  to  all  in  order  to  simplify  the  complex 
calculations.  In  actual  practice,  these  assumptions  do  have  their  own  effect  on  final  outcome.  One  of  such  assumption  is 
that  the  effective  coefficients  have  been  calculated  with  respect  to  average  electric  field  in  the  homogenized  piezoelectric 
composite.  In  actual  practice,  there  is  a  huge  difference  in  electric  field  in  fibre  and  matrix  phases  due  to  striking  contrast 
in  their  dielectric  properties.  Authors  [7-10]  have  used  these  same  set  of  assumptions  to  predict  overall  effective  properties 
of  PFRC.  Few  improvements  of  properties  have  also  been  reported,  by  doing  certain  modification  to  earlier  models.  Kumar 
and  Chakraborty  [9]  have  extended  study  of  Mallik  and  Ray  to  determine  the  overall  thermo-electro-elastic  properties, 
which  is  based  on  strength  of  materials  approach. 

Present  work  aims  at  developing  a  micromechanics-based  model  for  the  evaluation  of  effective  electro¬ 
mechanical  properties  of  PFRCs,  which  could  be  used  in  accurate  analysis  of  PFRCs  as  actuators  in  smart  applications. 
Micromechanics  model  for  predicting  the  effective  properties  of  the  piezoelectric  composite  is  based  upon  modified 
strength  of  materials  and  energy  method.  Results  obtained  has  been  compared  to  the  effective  properties  obtained  through 
strength  of  materials  method  [7-10]. 

2.  MICROMECHANICS  MODEL 

Piezoelectric  fibre  reinforced  ceramics  consists  of  piezoelectric  fibers  aligned  parallel  to  each  other  in  longitudinal 
direction,  surrounded  by  matrix  phase.  Bulk  materials  can  be  visualized  as  an  assembly  of  a  rectangular  representative 
volume  element  (RVE)  containing  both  fiber  and  matrix  which  can  be  repeated  in  three  dimensions  to  obtain  the  bulk 
material.  Figure  2  shows  the  transverse  cross-section  of  such  a  representative  volume  element  in  a  PFRC.  This  section 
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contains  micromechanics  modelling  of  such  RVE  with  modified  strength  of  material  method.  The  effective  coefficients 
obtained  using  modified  strength  element  would  be  same  as  that  of  RVE  in  average  sense. 


For  simplification,  few  assumptions  have  been  made  in  this  study.  It  is  assumed  here  that  the  composite  is 
subjected  to  a  constant  electric  field  in  the  direction  transverse  to  that  of  fibers,  electric  field  is  same  in  both  matrix  and 
fiber.  Piezoelectric  composite  is  uniform  throughout,  behaves  in  linear  elastic  manner,  and  there  is  no  slippage  between 
fibers  and  matrix.  Also,  fibers  are  continuous  and  aligned  parallel  in  longitudinal  direction.  Figure  1  shows  a  square  array 
of  fibers  of  circular  cross-section  and  figure  2  gives  representative  volume  element  (RVE)  of  such  an  array.  The  circular 
cross-section  of  fiber  is  converted  into  square  with  equivalent  area,  and  the  RVE  is  divided  into  sub  regions  a  shown  in 
figure  2.  Area  of  cross-section  of  original  fiber  is  a.f,  which  can  be  modified  to  an  equivalent  square  fiber  of  side  Sf.  Since 
length  of  side  of  RVE  is  s,  the  cross-section  area  of  RVE  becomes  s2.  As  cross-section  area  of  circular  fiber  is  equal  to 
cross-section  area  equivalent  square  fiber  i.e.  sf  =  dj-  or  .sy  =  ^faf-  Fiber  volume  fraction  for  longitudinal  fiber  reinforced 
composite  is  the  ratio  of  cross-section  area  of  the  fiber  to  cross-section  area  of  RVE  i.e. 


(1) 


(2) 


As  observed  from  figure  2  matrix  dimension  is  sm  =  s  —  Sf,  from  eqn.  (2)  it  can  be  written  as 


(3) 


z,3 


Figure  1:  Schematic  Diagram  of  PFRC  Lamina. 


Figure  2  illustrates  that  RVE  can  be  subdivided  into  sub  regions  I  and  II,  region  I  is  matrix  dominated  while  II 
is  a  composite  consisting  both  fiber  and  matrix  phases.  In  this  approach,  firstly  effective  properties  of  sub  region  II  are 
computed  as  series  connection  of  fiber  and  matrix  subjected  to  transverse  normal  stress,  which  is  followed  by 
calculation  of  overall  effective  properties  taking  sub  region  I  and  II  as  parallel  connection  subjected  to  transverse 
normal  stress. 
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Constitutive  equations  for  transversely  isotropic  piezoelectric  active  fibers  in  Voigt’s  two  index  notation  are  [7-10]. 
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The  constitutive  equations  for  isotropic  passive  matrix  phases  in  Voigt’s  two  index  notation  are 
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here  C[r  e*-,  are  elastic  constants,  piezoelectric  constants  and  dielectric  constants,  respectively  and  ajj,  e£h 
E£,  Df  are  the  stress,  strain,  electric  field,  electric  displacement,  respectively.  Superscript  x  denotes  piezoelectric  fiber  (f) 
or  matrix  (m)  phases,  respectively. 


The  resultant  constitutive  equations  of  composite  can  be  written  as 
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(7) 


Where, Cy,  efj,  tcfj  are  to  be  determined  in  order  to  evaluate  complete  electro-mechanical  behaviour  of  the 
composite. 


Considering  series  connection  of  fiber  and  matrix  in  sub  region  II  subjected  to  transverse  loading,  geometric 
compatibility  requires 


8&=8"I  +  SiT 


(8) 


S22  ,  Sir  are  displacements  in  fibers  and  matrix  when  transverse  loading  in  direction-2  is  applied.  Consequently, 


*22$  =  42sf  +  £22Sm 

Using  Eqns.  (1),  (2)  Eqn.  (9)  can  be  written  as 
*22  =  sjVf42  +  (1  “  44^22 

Similarly,  for  direction-3,  the  relation  can  be  written  as 

*33  =  ^4 3  +  c1  - 

As  perfect  bonding  is  assumed,  therefore  the  strains  are  same  in  direction-2 

e[\  =  e{j_  =  e'T  (12) 

Also,  an  electric  field  is  applied  in  transverse  direction  only  and  is  same  in  both  fiber  and  matrix  for  sub  region  II 
11  —  rf  — 


(9) 


(10) 


(11) 


p“  _  pJ  _  pm 
c3  ~  c3  —  C3 


(13) 


Stresses  in  fiber  and  matrix  phases  in  transverse  directions  2  and  3  are  equal  and  written  as 

(14) 

33  ~  u33  (15) 

Using  Eqns.  (4),  (5),  (10),  (11),  (12),  and  (13)  into  Eqn.  (14)  and  (15)  the  following  Relations  are  obtained 

Axe[2  +  A2 4  =  Bre{x  +  B2ec22  +  B3e$3  -  B4E§  (16) 


4l  =  °22 
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www.tjprc.org 


SCOPUS  Indexed  Journal 


editor@tjprc.  org 


122  Sanjeev  Singh  &  Saroja  Kanta  Panda 

A2ef22  +  A3e'3  =  Bsecu  +  B3ec22  +  B6ec33  -  B7EC3  (17) 
where, 

Ai  =  (l~  JvfYL  +  JvfCPi 
a2  =  (!  -  -Jvf)cL  +  4vfC£ 
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On  solving  Eqns.  (16)  and  (17),  the  following  Eqns.  are  obtained 

=  fell  +  22  +  heTi  ~  UE3  (19) 

eL  =  9ieii  +  92^22  +  93^33  ~  94E3  (20) 

where. 
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Using  Eqns.  (4),  (5),  (11),  (12),  (13)  into  Eqns.  (19)  and  (20)  transverse  stresses  in  2-direction  can  be  expressed  as 
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a22  ~  +  ^2e22  +  h£33  U^3  (22) 

+  m2ec22  +  m3e^3  -  m4E%  (23) 

Where, 
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Subsequently,  considering  sub  regions  I  and  II  as  combined  unit  in  RVE  and  applying  following  conditions 

°2C2  =  +  (!  -  22  (24) 
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The  effective  stiffness  and  piezoelectric  coefficients  obtained  as 
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e33  —  e33  ^23/4  +  ^33^4 

3.  ENERGY  METHOD 

Energy  functions  for  piezoelectric  composite,  fiber  and  matrix  are  given  by 


uc  =  J  | -i(ec7CcV  +  2£CW )  + 

(26) 

Uf  =  /{-i(V  W  +  2^W)  +  1-EfTefTEf}dV 

(27) 

Um  =  f\-^(emTCmTern  +  2EmTemTem)  +  ^EmTemTEm}dV 

(28) 

Under  the  given  state  of  stress,  total  energy  stored  in  composite  is 

matrix.  Thus, 

equal  to  sum  of  energies  stored  in  fibers  and 

uc  =  uf  +  Um 

(29) 

Stresses  and  strains  in  fibers  and  matrix  are  defined  in  terms  of  corresponding  composite  quantities,  as  given 

below 
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Ef  =  gwEcz  (47) 
Using  Eqns.  (30)-(35),  (46)  and  (47)  following  relations  are  obtained 

eii  =  «feici  +  a2e22  +  4e33  -  a{E%  (48) 

e?i  =  c1±  +  a?ec22  +  a™ec33  -  <£3C  (49) 
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Using  Eqns.  (30)-(49)  into  Eqn.  (29),  equating  the  coefficients  of  like  terms  both  side,  effective  coefficients  are 

obtained. 


rc  ■ 

rf  rm 
_  LnLli 

9icii+92c[i 

rc  - 

=  a{  9?,  c{2  +  a?g4C£ 

rc  ■ 

=  a{gsc[3  +  off  g6C^ 

rc 
u22  ‘ 

=  ■  f  ^  ^  cm  ^  +  gic22  +  9lc?2  +  2^53  <7/2  + 

C11  L11 

rc 

u23  ‘ 

=  af3S3E(2  +  a™g4C%  +  af2g5c{3  +  a?g6C£  +  g3g5C^3  +  ,g4,g6C2^ 

E33 ' 

(giCi2-g5c{2)  ( 92Ci2-96c{2 )  2rf  1  2rm  ,  n  /  rf  \  n  m  rm 

=  f  +  rm  +  55C33  +  9e  +  2  aJ3gsC(3  +  2a^g6C^ 

e31  : 

=  a^i7e{i  +  <5i8e™i 
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e3i  =  a2  517^31  +  arises!  +  939tie{2  +  g^wefz  +  a{g3c{2  +  a^g4C^2 

el3  =  a3  $i7e;fi  +  +  9s9nel3  +  ^  +  af4gsc[3  +  a^g6C^3 

It  should  be  noted  that  constants  gt  (t  =  1,2,3 . 18)  are  experimental  parameters  and  must  be  obtained  through 

experimental  techniques.  However,  experimental  determination  of  these  parameters  is  out  of  the  scope  for  the  present 
study.  Thus,  these  parameters  are  obtained  with  the  theoretical  assumptions  of  modified  strength  of  materials  approach  as 
discussed  in  the  previous  section. 

4.  RESULTS  AND  DISCUSSIONS 

Above  derived  equations  obtained  through  modified  strength  of  materials  and  energy  methods  has  been  used  to  study  the 

diversion  of  effective  properties  of  PFRC,  obtained  by  strength  of  materials  method  [7-10].  Table  1-2  shows  the  elastic, 

electric  and  dielectric  properties  of  the  constituents  i.e  fiber  and  matrix  phases. 


Table  1:  Elastic  Properties  of  Piezoelectric  Fibers  and  Matrix  [9,10] 


Material 

Cn(GPa) 

C12(GPa) 

Ci,(GPa) 

C  33(GPa) 

C44(GPa) 

Cfifi(GPa) 

PZT-5 

121 

75.4 

75.2 

111 

22.8 

21.1 

Epoxy 

3.86 

2.57 

2.57 

3.86 

0.64 

0.64 

Table  2:  Piezoelectric  and  Dielectric  Properties  of  Fibers  and  Matrix[9-10] 


Material 

e3i(C/m2) 

e33(C/m2) 

e15(C/m2) 

/c11(109C/Vm) 

k33(10  9  C/Vm) 

PZT-5 

-5.4 

15.8 

12.3 

8.11 

7.35 

Epoxy 

0 

0 

0 

0.0797 

0.0797 

The  study  carried  here  provides  a  comparative  study  of  variation  in  effective  properties  as  a  function  of  fiber 
volume  fraction  and  matrix  properties  among  strength  of  materials,  modified  strength  of  materials  and  energy  method.  The 
following  non-dimensional  parameters  have  been  used  in  the  analysis. 


R 


32 


e 

e 


c 

32 

/ 

32 


R33 


e 

e 


c 

33 

/ 

33 


Figures  3-6  show  the  comparative  study  of  effective  elastic  coefficient  of  composite  with  respect  to  change  in 
fiber  volume  fraction.  Results  clearly  shows  that  modified  strength  of  method  and  energy  method  are  close  in  values,  while 
strength  of  materials  overestimates  the  effective  coefficient  Cn  while  it  underestimates  effective  coefficient  C12.  Figure  5 
shows  effective  coefficient  estimated  from  strength  of  materials  method  is  in  close  proximity  with  modified  strength 
method  and  energy  method.  Again,  a  huge  underestimation  of  effective  coefficient  C33  is  reported  as  shown  in  figure  6. 
These  wide  deviations  of  effective  properties  occur  due  to  fact  that  in  strength  of  material  approach,  there  has  been  certain 
assumptions  which  provide  very  poor  results.  In  strength  of  materials  approach  [7-10]  it  is  assumed  that,  equal  strains  are 
produced  in  fiber  as  well  as  matrix  phases  when  stress  is  applied,  while  in  actual  practice  these  strains  ought  to  be 
different.  Modified  strength  of  materials  approach  refines  result  for  transverse  stresses,  while  energy  approach  also  gives 
better  result,  as  Poisson’s  ratio  mismatch  is  tackled  to  a  great  degree  through  this  model. 
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Figure  3:  Variation  of  Elastic  Coefficient  Cn  of  PZT-5/Epoxy 
Composite  with  Fiber  Volume  Fraction. 


Figure  4:  Variation  of  Elastic  Coefficient  C12  of  PZT-5/Epoxy 
Composite  with  Fiber  Volume  Fraction. 


Figure  5:  Variation  of  Elastic  Coefficient  C13  of  PZT-5/Epoxy 
Composite  with  Fiber  Volume  Fraction. 
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Figure  6:  Variation  of  Elastic  Coefficient  C33  of  PZT-5/Epoxy 
Composite  with  Fiber  Volume  Fraction. 

Similar  trends  of  overestimation  and  underestimation  is  also  reported  for  R31,  R32,  and  R33  for  strength  of 
materials  in  comparison  with  modified  strength  of  materials  and  energy  approach.  Figure  7  shows  value  of  ratio  of 
piezoelectric  coefficient  R31  which  has  been  highly  overestimated  in  strength  of  materials  method  [7-10]  as  compared  with 
modified  strength  of  materials  and  energy  method.  While  results  shown  in  figure  8  Indicates  the  values  of  R32  are 
somewhat  in  good  agreement  with  each  other,  at  extremely  lower  and  higher  values  of  fiber  volume  fraction.  In  mid  ranges 
of  fiber  volume  fraction,  there  is  symmetric  divergence  of  properties  between  these  three  methods.  In  this  case,  strength  of 
materials  approach  and  modified  strength  of  materials  approach  forms  lower  and  upper  bound  of  energy  approach, 
respectively.  Figure  9  shows  ratio  R33  is  in  good  agreement  for  strength  of  materials  method  and  energy  method,  while 
there  is  a  divergence  of  modified  strength  of  materials  results  from  these  two  methods. 


f 

Figure  7:  Variation  of  Ratio  of  Piezoelectric  Coefficient  -  e3l/ e31 
PZT-5/Epoxy  Composite  with  Fiber  Volume  Fraction. 
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f 

Figure  8:  Variation  of  Ratio  of  Piezoelectric  Coefficient  e32/e32  PZT- 
5/Epoxy  Composite  with  Fiber  Volume  Fraction. 


f 

Figure  9:  Variation  of  Ratio  of  Piezoelectric  Coefficient  e33/ e33  PZT- 
5/Epoxy  Composite  with  Fiber  Volume  Fraction. 

Though  there  is  a  huge  variation  in  dielectric  properties  of  the  constituent  phases  of  PFRC,  figure  10  shows  that 
the  overall  dielectric  coefficient  of  composites  is  in  good  agreement  with  each  other  for  all  three  methods  i.e.  strength  of 
materials,  modified  strength  of  materials  and  energy  method. 


Figure  10:  Variation  of  Dielectric  Coefficient  k33  PZT- 
5/Epoxy  Composite  with  Fiber  Volume  Fraction. 
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5.  CONCLUSIONS 

In  this  paper,  two  distinguished  method  of  micromechanics  i.e.  modified  strength  of  materials  and  energy  method  have 
been  used  to  study  the  effective  properties  of  piezo  fiber  reinforced  composites.  The  results  have  been  used  to  analyse 
deviation  of  effective  properties  estimated  in  earlier  works  [7-10]  using  strength  of  materials  method.  In  various  works 
[7-10],  the  strength  of  materials  method  models  the  problem  on  assumptions  that  the  equal  strains  are  produced  in  both 
fiber  and  matrix  phases,  while  this  assumption  doesn’t  have  significant  effects  on  results  while  the  longitudinal  stress  is 
applied  to  the  composite.  In  case  of  transverse  loading,  the  results  obtained  from  strength  of  material  would  be  far  from 
accurate.  Modified  strength  of  materials  method  refines  the  result  obtained  from  strength  of  materials  method  and  provides 
a  better  estimated  effective  property.  Energy  method  modelling  of  the  problem  eliminated  the  effect  of  Poisson’s  ratio 
mismatch  between  fiber  and  matrix  phases.  Thus,  present  study  reveals  better  models  for  micromechanics  analysis  than  the 
reported  strength  of  materials  methods  for  various  coefficients,  and  expected  to  give  more  precise  results  than  existing 
method. 


It  makes  tailoring  of  properties  of  piezoelectric  composites  that  contains  piezoelectric  fibers  in  piezoelectrically- 
inactive  matrix,  having  superiors  and  improved  effective  electrical  and  mechanical  properties.  In  addition,  these  models 
could  be  served  as  a  benchmark  for  numerical  analysis  of  PFRC  problems  and  can  be  utilised  to  model  coupled  more 
piezoelectric  composite  problems,  subjected  to  more  complex  boundary  conditions. 
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